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Abstract 

A combined fit is presented to data on pp annihilation in flight to final states tjt/^ti^, TiPn^, rp], rp]' and n'n'^. The 
emphasis lies in improving an earlier study of rjif'Tf' by fitting data at nine p momenta simultaneously and with 
parameters consistent with the two-body channels. There is evidence for all of the / = 0, C = +1 states expected 
in this mass range. New resonances are reported with masses and widths {M,T) as follows: Z^'- = (2328 ± 38, 
r = 240 + 90) MeV, 1++ (1971 + 15, 240 + 45) MeV, 0"+ (2285+20, 325 + 30) MeV, and 0"+ (2010+^,f,, 270 ±60) MeV. 
Errors on the masses and widths of other resonances are also reduced substantially. All states lie close to parallel 
straight line trajectories of excitation number v. mass squared. 



Data from pp interactions in flight have the great merit of allowing a direct study of ^-channel meson resonances 
R in formation reactions of the type pp ^ R ^ A + B. Many decay channels A, B may be studied. Earlier, we 
have presented data on pp rjTfiT^ [1], in which evidence was found for a number of / = C = +1 mesons. 
That analysis fitted partial waves separately at each of nine p momenta; resulting magnitudes and phases were then 
interpreted in terms of resonances. [Table 1 shows the relation between beam momenta and centre of mass energies.] 
Here our objective is to fit these data at all momenta simultaneously. The fit to resonance parameters is also made 
simultaneously to data on pp jfin'^, rjTj, rjTj' and tt^tt^ reported earlier [2,3]. 

This combined fit is constrained strongly by the requirement of consistency between channels and gives accurate 
determinations of masses and widths of most resonances found in the data. Fluctuations which were present in 
Refs. [1] at individual momenta are eliminated; this is particularly important in restricting phases to a smooth energy 
dependent behaviour consistent with analyticity. In the earlier work, fluctuations in the large high partial waves 
obscured details in small ampUtudes, particularly low partial waves. We are now able to identify further ^-channel 
resonances with 7^'' = 4"^, 1^^ and 0"^. The essential results of the new analysis are summarised by the masses 
and widths given in Table 2. These results supercede earlier determinations in Refs. [1] and [3]. The present paper 
completes our analyses of channels with 7 = 0, C = +1. 

We begin by outlining the considerations which influence the form of our partial wave analysis. The first point 
is that the well-known /4(2050) appears strongly in all data sets, although we find a distinctly lower mass than that 
quoted by the Particle Data Group (PDG) [4]. The /4(2050) acts as an interferomenter. Its interferences determine 
relative phases of other partial waves. These relative phases are found to vary little with mass, hence requiring that all 
partial waves follow a similar resonant behaviour to /4(2050). The way the rjnn analysis goes is that one immediately 
finds strong J^^ = 3"^"^ and 2""^ resonances. The two-body data require 2^^ and O"^"^ resonances. Then, with increasing 
attention to detail, 4"^, 1^^ and 0"^ resonances emerge from the rjnn data. 

We therefore parametrise each partial wave amplitude as a sum of resonances, plus a constant or slowly varying 
background where necessary. In analysing earlier data [1,3], we have found that resonances mostly cluster into two 
groups (a) from 1920 to 2050 MeV around /4(2050), (b) from 2220 to 2320 MeV, around /4(2300). 

The background terms may parametrise the tails of resonances below the pp threshold or the eff'ects of f-channel 
exchanges. Singularities due to those exchanges are distant {s <\ GeV^), so it is to be expected that backgroimds wiU 
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Beam momentum 
(MeV/c) 



(MeV) 



600 
900 
1050 
1200 
1350 
1525 
1642 
1800 
1940 



1962 
2049 
2098 
2149 
2201 
2263 
2304 
2360 
2409 



Table 1: Momenta at which data are available and corresponding centre of mass energies. 



contribute mostly to low partial waves. This is what we find. The strong high partial waves with J^^ = 4"^'^ and 3"^"^ 
are consistent with no background. In low partial waves (J < 2), the background is parametrised as a broad resonance 
(in most cases below the pp threshold) or as a constant. Parametrising with resonances or constants guarantees that 
partial wave amplitudes obey the important constraint of analyticity, since a Breit-Wigner amplitude is analytic. We 
see no evidence that strong threshold efifects are present to perturb such a parametrisation. 
Each partial wave amplitude then takes the form: 



The factors B{pp) and B{AB) are standard Blatt-Weisskopf centrifugal barrier factors which guarantee the correct 
threshold behaviour for each channel; expressions are given in Ref. [5]. A common radius for the centrifugal barrier 
is fitted to all partial waves. There is a strong optimum at 0.829 ± 0.021 fm. The full widths F, of all resonances 
are taken to be constant because of the large number of open channels. Each resonance is fitted with a real coupling 
constant gi and a phase 0, . 

We now discuss the channels fitted to the rjnn data. Fig. 1 shows mass projections on to nn and rfji at two 
representative momenta; further figures are to be found in Refs. [1]. In Fig. 1, data are uncorrected for (small) 
variations of acceptance, which are included in the maximum likelihood fit. Histograms show the results of the 
present fit. In Ref [1], final states fitted to rinPif' data were dominantly fl2(1320)7r, /2(1270)77 and rjcr, where cr stands 
for the nn S-wave amplitude. There were small additional contributions from ao(980);r, /o(980)77 and /o(1500)77. In 
the rjn mass projection at the higher beam momenta there is a distinct shoulder in the mass range around 1450 MeV. A 
fit to this part of the the mass spectrum and higher rjn masses requires further small contributions from ao(1450)7r and 
a2(1660);r. The evidence for a2(1660) Jjn from the Crystal Barrel experiment is given in Ref. [6]. Parameters for 
ao(1450) are fixed to values of the PDG. However, evidence presented here for s-channel resonances does not depend 
significantly on these small a2(1660)7r or ao(1450) contributions. We have tried including ^)(1370) in the fit, but find 
neghgible eff'ect. It is known to decay weakly to nn [7] compared with /o(1500), and the latter is already a small efi'ect 
in the present data. 

The 2-body channels nn, rjrj and rjrj' are related by SU(3), since n, rj and are members of a single nonet. 
Formulae for these constraints are given in Refs. [3]. In outhne, SU(3) allows 7777 and 7777' amphtudes to be predicted 
from nn using the well-known composition of 77 and 77' in terms of {uu + dd)/ "V^ and ss. Here we use the assumption 
that pp does not couple directly to ss; this assumption has been tested in Refs. [3] and is well obeyed with one 
striking exception. The /o(2105) decays more strongly to 7777 than to nn by a factor 1.88; this compares with the SU(3) 
prediction of 0.8"^ = 0.42. That is, the decay to 7777 is a factor 4.5 stronger than predicted by the SU(3) relation. Since 
the /()(2105) is produced strongly, this feature suggests it has exotic character. Its coupling constants to nn, 7777 and 7777' 
are therefore fitted freely. 

For 77;r;r, amphtudes for a2(1320);r and /2(1270)77 are again in principle related by SU(2) constraints if one assumes 
ideal mixing for the 2"^ nonet. These have been tried in the fit, but do not work well. This is probably because (a) 
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there is a large mass difference between these channels and (b) the f2(l270)T] threshold is close to the region we are 

analysing. Consequently, momenta in these channels are substantially different. This is likely to have large effects on 
matrix elements. We find it necessary to fit magnitudes and phases of a2il320)n and /2(1270)?7 amplitudes freely. 

Mass differences between n, tj and tj' also give rise to differences in momenta q in nn, rjrj and tjtj' channels, hence 
significantly different form factors and centrifugal barrier effects. We fit a form factor exp(-Qr^^) to all decays. The 
value of a optimises at 0.92 + 0.10 GeV"^. 

We find that two minor improvements may be made to the fit. The phases 0, of each resonance arise from multiple 
scattering in initial and final states, as explained in Ref. [8]. Firstly, to allow for overlap of resonances (hence departure 
from strict Breit-Wigner forms), phases <pi are allowed to vary by up to ±20° for 7777 and 7777' with respect to the Tin 
value. Secondly, we find significant improvements by allowing coupling constants g, to vary for 7777 and 7777' from 
their SU(3) values by factors constrained to the range 0.7 to 1.3. This probably reflects differences in form factors 
for different matrix elements, which may be affected by masses, hence momenta, in the two-body final state. These 
refinements have only minor effects; without them, conclusions on fitted resonance masses and widths change by only 
a few MeV. 

For = 2^ and 4^, pp may couple with orbital angular momentum L = 7 ± 1 (e.g. and ^Fa)- An ideal reso- 
nance should have the same phase for both L values (via multiple scattering through the resonance to each channel). 
We therefore take the ratio of coupling constants rj - gj+\ lgj-\ to be real. In fitting rjjui data, where a final state such 
as /2( 1270)77 may also have two or more L values, the ratio of coupling constant is likewise constrained to be real 
in initial fits. Again, we find small but significant improvements to the fit if these ratios for /277 and a27i decays are 
allowed to depart from real values by phase angles in the range up to ±20°. But there is little effect on fitted masses 
and widths, merely a small improvement in the quality of the fit to data. 

In Refs. [3], data on mi, rjrj and 7777' required four 7^'' - 2^^ resonances. In addition, data on pp — > {rjTffjf' 
provide evidence for a broad 2^ 7777 contribution with M - 1980 ± 50, Y - 500 ±100 MeV [9]. Central production 
of An provides even clearer evidence with very similar mass and width [10]. We find that this broad component 
improves the present combined fit to ?77r;7r and 2-body channels strongly. Here we find a clear optimum for the mass at 
M = 2010 + 40MeVandforfliewidfliat495 + 50MeV. These values are consistent with Refs. [9] and [10]. However, 
those earlier determinations have the advantage that the broad 2^ component is clearly visible by eye, hence providing 
the incentive for trying it here. Without this component, large interferences develop between the four 2*' states, so as 
to simulate this broad contribution. 

We have searched for ambiguous solutions by (i) removing each resonance (or background) one by one and re- 
optimising the rest, (ii) changing signs of ampUtude ratios r one by one and (iii) moving resonances in steps of 20 
MeV in mass and 40 MeV in width over large ranges (9 steps each). Each iteration of the overall fit takes 10 minutes 
of computing and one solution converges in typically 20-50 iterations. So ~ 500 altematives have been explored. All 
variants collapse back to the same solution. 

Intensities of all amplitudes fitted to t^t^tiP are displayed in Fig. 2. They have been corrected for all 77 and 71^ de- 
cays, i.e. for the 39.25% branching fraction of 77 —> yy and the 98.798% branching fraction of 7^ yy. Typical errors 
are ±5% - 10% for the larger partial waves, increasing to +35% for the smallest. Amplitudes which contribute less 
than 0.3% of the integrated cross section are found to change log likelihood by < 40 and are omitted.Table 2 shows 
masses and widths of fitted resonances. Statistical errors are neghgible, so quoted errors cover the range of solutions 
observed when components of the fit are varied. It is apparent, for example, that there are small systematic discrepan- 
cies between data on final states n'n'^ and n^n^, and some masses and widths change by a few MeV according to the 
way one weights different data sets. Errors include systematic changes when the radius of the centrifugal barrier is 
varied within its error The last column of the Table shows changes in log likelihood when each component is removed 
from the fit to rjn^n^ data and others are re-optimised. Our definition of log likelihood is such that it changes by 0.5 
when one parameter changes by one standard deviation. For 2"^ and 4"^ components, there are also large changes in 

arising in the fit to two-body data; those changes are close to values quoted in Refs. [3]. Consequently, all of the 
components in the fit are highly significant. The last 6 lines of Table 1 show, for completeness, parameters of / = 1 
resonance fitted to two-body data; there are small changes from values of Refs. |3J. 

The fit to nn, rfrj and 7777' has changed little from that of Refs. [3], so attention will be concentrated here on fits 
to rjnn. For those data there are some significant changes from the fit described in Refs. [1]. The main origin of 
these changes is that ratios r of amplitudes with L = 7 ± 1 are defined well by the two-body data, specifically by 
the polarisation data on pp — > ;r~;r"^. Using these ratios in the combined fit introduces significant changes to 2*' 
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Table 2: Resonances fitted to the data. Errors cover systematic variations observed in a variety of fits with difl'erent ingredients. Values in 
parentheses are fixed from other data [18]. The sixth colunm shows changes in S = log likelihood when each component is omitted in turn from 
the fit to rjifiTfi data and remaining components are re-optimised. The final column highlights those resonances which are new. 
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amplitudes in rjnn. A knock-on effect is that there is also a large improvement in the fitted 3"^ amplitude towards the 
top of the mass range. 

Argand diagrams are displayed in Fig. 3. Crosses mark individual beam momenta; Table 1 lists these beam 
momenta and the corresponding centre of mass energies. In many cases, one discerns in Fig. 3 that the amplitude 
varies rapidly in the low mass region and again in the high mass region. This indicates the presence of two resonances. 
However, only in Figs. 3(h) and (e) does one see conspicuously separate loops due to different resonances. 

The data require a larger and more Unear phase variation than can be produced by a single resonance. A sequence 
of resonances conspires to produce approximately a linear phase variation with mass squared s. In electronics, it 
is well known that a linear phase variation with frequency may be obtained with a Bessel filter [11], which has a 
sequence of poles almost linearly spaced in frequency. Meson resonances appear to behave in an analogous way as a 
function of i = M^. Physically, such a system generates a time delay which is independent of frequency, since group 
velocity is proportional to the gradient of phase v. frequency. 

We shall now comment on each 7^ in turn, beginning with the high partial waves which are most conspicuous. 
For 4^, the different peaks in Figs. 2(g) and (h) for /2(1270)7/ (full curves) and fl2(1320)7r (dashed) obviously demand 
two resonances. The lower one optimises at a mass of 2018 ± 6 MeV, considerably lower than the PDG average of 
2044 ±11 MeV. The effect of the centrifugal barrier is large and the intensity of the ^Fa, partial wave peaks at 2080 
MeV in Fig. 2(g). Differences in earlier determinations of the mass may well depend on varying treatments of the 
centrifugal barrier. Our mass is determined essentially by the maximum in the 'speed plot', i.e. in the movement of 
the amplitude in the Argand diagram; it is found to be quite insensitive to the radius chosen for the centrifugal barrier. 
As the radius increases, all masses go up; however, with our error for the radius of the barrier, the contribution to the 
uncertainty in mass is only 0.7 MeV. The width of the resonance does, however, correlate much more strongly with 
the radius, and this contributes an uncertainty to the width of 3.3 MeV. The L = 3 centrifugal barrier in the pp channel 
suppresses strongly the coupling of the lower 4^ resonance to /2(1270)/;. The lower resonance couples purely to pp 
^F4, while the upper one requires a large component. 

For 7^ = 3^, the different peaks in Fig. 2(d) for f2r] (full curves) and (dashed) again clearly demand two reso- 
nances. The Argand diagram of Fig. 3(h) also clearly demands two resonances. The upper resonance is considerably 
stronger in the present fit than in Refs. [1]. 

For j'' = 2"^, the rjnn data alone do not resolve /2(1934) and /2(2001) clearly, because both lie at the bottom end 
of the available mass range. One sees from Table 2 that removing the /2(2001) produces a change in log hkehhood in 
T]nn data alone of only 168. This is because it may be interchanged to some extent with contributions from /2(1934). 
For /2(1934) and /2(2001), the two-body data play a very important role, for two reasons. Firstly there are extensive 
data on the tt tt^ channel down to 360 MeV/c (a mass of 1910 MeV), which determine quite well the mass and width 
of the /2(1934). This state is consistent in mass with the the /2(1920) of GAMS [12] and VES [13], though they find a 
narrower width. Secondly, the polarisation data provide a sensitive determination of the ratios of amphtudes between 

and ^P2. The values of r2 given in Table 2 are for ratios of coupling constants ^F2/^P2, i-e. after factoring out the 
effect of the centrifugal barriers. Both /2(1934) and /2(2001) are definitely required by the two-body data. If either of 
them is removed from the fit,^^ for 2-body data increases by > 2000, as reported in Refs. [3]. 

The earlier analysis of rjnn data [1] showed the requirement for three 2"^ states at 2020 ± 50, 2240 ± 40 and 
2370 ± 50 MeV. These masses now adjust naturally by small amounts to those of Table 2. The lowest 2"^ state at 1934 
MeV contributes strongly to pp — > nn and is dominantly 2; that at 2001 MeV is largely ^^2 and is close in mass 
to and resonances. The /2(2240) and /2(2293) are dominantly ^^2 and ^^"2 respectively. Each of the Argand 
diagrams of Figs. 3(k) lo (o) show at most a requirement for two 2"^ states. However, the two-body data require four 
2"^ states, as reported in Refs. [3]. The /2(1934) lies below the available mass range for rjnn data. 

We have speculated earUer on the possibility that the broad /2(1980) has a large component due to the 2+ glueball 
[14,15]. In this context, its relative coupling to nn and 7777 is important. We fit these freely, with the result g^^/g^^o - 
0.72 ± 0.06. This compares with the prediction 0.41 for (uii + dd)/ and 1 for a glueball. The result Ues midway 
between the two. In Ref. [3], the amplitude was expressed in terms of flavour mixing to a linear combination of states 
cosOImm + dd > I yll + sinOIss >. With this parametrisation, the flavour mixing angle O optimises at 23.6 ± 3.5° 
compared with the value 35.6° for a glueball and 0° for {uu + dd). Some mixing between a glueball and neighbouring 
qq states is likely. 

For 7^ = 0^, there are only minor changes from Ref. [3]. The /o(2105) makes a large and very well defined 
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contribution to t]T], and a small contribution to nn. It requires a mixing angle O = (58 ± 5)°. Its strong production 

and large flavour mixing angle suggests exotic character. It is clearly an 'extra', non-<7^ state. Its mass optimises at 
2102 + 13 MeV and its width at 21 1 +29 MeV. These values agree closely with those found in the E760 experiment 
[16]. Since their determinations of the width is somewhat more precise, namely 203 ±10 MeV, we use this value in 
the final fit. The data also require the presence of another nearby broad /o(2040), in good agreement with WA102 
parameters [17] and a further /o(2337) at higher mass. 

For J'''^ = 2""^, earUer data on pp — > ■qT^iifn'^ demonstrate the presence of two states at 1860 and 2030 MeV 
[18]. The first of these decays largely to [/277]l=o and the second decays weakly to this channel. The present Tjnn data 
extend only down to 1960 MeV and therefore do not resolve these two states. We therefore fix their masses, widths 
and branching ratios between /277, and UQn channels to the values of Ref. [18]. The [/2?7]l=o channel makes a 
dominant contribution to the 2"^ rjnn partial waves at low mass, see Fig. 2(c), full curve. The presence of a further 
772 at high masses is seen most clearly in 'D2 — > [f2^]L=2, Fig. 2(c) chain curve; it is also visible in the smaller 
components [ao(980);r]z,=2, Fig. 2(k) dashed curve and in 'D2 — > [/o(1500)77]i=2, Fig. 2(k) chain curve. Its mass, 
2267 ± 14 MeV, is somewhat lower than the earlier determination [1], 2300 ± 40 MeV. However, it agrees closely 
with a conspicuous peak observed [19] in the integrated cross section for pp — > ti't^'tiP. There, a mass of 2248 + 20 
MeV is found. The width from present data is 290 + 50 MeV. However, the width is determined better by rj'n'^n^ data: 
280 ± 20 MeV, and we fix it at this value in the final fit to present data. 

We now come to the new states. There is a well known ps (qq ^Gs) listed by the PDG at 2330 MeV. One expects 
a qq 4"^ state (qq 'G4) of similar mass. We now observe small but well determined contributions from a 4"^ state 
at 2328 + 38 MeV in 'G4 — > [fl2^]L=2^ Io2^]l=4 and [aon]i=4, Fig. 2(m). This resonance, although it contributes the 
smallest change in log likelihood in Table 2, is very stable throughout all fits. There is no doubt of its presence, despite 
the fact that it contributes only 2% of the integrated cross section. However, there is quite a large error on the width. 

Next we consider 1^^. The earlier work of Ref. [1] found a strong [a27i]L=i intensity near the pp threshold. In the 
new combined fit, the 1^^ amplitude is large in several channels around 1975 MeV, falling rapidly at higher masses, 
see Fig. 2(j). Despite the fact that it lies at the bottom end of the available mass range, parameters of the resonance 
in this range are very stable: M = 1971 ± 15 MeV, F = 241 + 45 MeV. Argand diagrams of Figs. 3(c) and (d) for 
[/2'?Jl=i and [o-ri]L=\ require rapidly varying phases at low mass when one remembers that the amphtude goes to zero 
at the pp threshold; this rapid phase variation is a clear signature of the resonance. 

At high masses, the phase variation observed in Refs. [1] suggested a l"^"^ resonance at 2340 + 40 MeV with 
r = 340 ± 40 MeV. We now find that this second l"^"^ resonance is definitely required. Without it, log likelihood is 
worse by 882; statistically this is a 25 standard deviation effect when one allows for the number of fitted parameters. 
Nonetheless, the mass of this state is the least well determined of all the resonances: M - 2310 ± 60 MeV. The 
reason is that it makes small contributions to the distinctive channels a2(1320);r and /2(1270)77, and there are large 
interferences with the 0" amplitude. As the resonance mass is changed, these interferences are able to absorb the 
variation with small changes in log UkeUhood. This resonance is most clearly visible in [/2?7]l=i, Fig. 2(b), full curve. 

Lastiy we consider 0""^. Here there is a large change, compared with Refs. [1], in the fit to the broad background 
in the Dalitz plot. Earlier this was fitted purely with the rjcr channel, although it was evident in Refs. [1] that the 
rjn mass spectrum was not fitted perfectly around a mass of 1450-1650 MeV. For that reason, conclusions were not 
drawn earlier about structure in the 0" partial wave. This feature has now been improved substantially by the addition 
of channels 00(1450)7?" and a2(1660)7r0. 

The at] and /o(1500)77 channels, Fig. 2(i) (fuU and dotted curves) now definitely require the presence of a high 
mass resonance with M = 2285 ± 20 MeV. In addition, there is a very strong contribution close to the pp threshold in 
several channels of Figs. 2(a) and (i). A resonance is definitely required in this mass region. Without it, log likelihood 
gets worse by a very large amount, 1189. The effect of dropping it is illustrated by the dashed curves of Figs. 3(a) 
and (b). For this partial wave, amphtudes do not go to zero at the pp threshold, but instead go to values described by 
scattering lengths. It is not possible to exclude the possibility that the pp threshold plays a strong role (a cusp at the 
pp threshold). This makes the identification of the mass of the lower resonance difficult. The optimum is at IQIQ"^^ 
MeV. 

There is just one significant change in the fit to two-body data, compared with Refs. [3|. There is a strong 3 
resonance in the low mass range. It interferes strongly with 2^* states in data on pp n^n^. Small changes in the 
masses and widths fitted to the 2"^ states have had the effect of increasing the mass of this pa from 1960 ± 15 MeV [3] 
to 1981 ± 14 MeV; the fitted width has also increased sUghtly. 
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Fig. 4 illustrates the quality of the fit at two momenta. It shows the angular distributions for production of mass 

regions centred on (a) and (d) /2(1270), (b) and (e) fl2(1320), (c) and (f) flo(980). In all cases, there are background 
contributions underneath these resonances, particularly in the case of the small ao(980) signal; these backgrounds are 
included in Fig. 4. The angular distributions are well reproduced by the fit in all cases. 

We now turn to the interpretation of the results. Fig. 5 shows plots of mass squared v. excitation for all . In 
making this plot, we use the K-matrix mass of 1598 MeV for /2(1565) [20], and 1400 MeV for /o(1370) [7]. The 
reason for this choice is that the well known linear mass relation between A(1230), S(1385), H(1520) and Q"(1672) 
works well when one uses for the A the mass at which the phase shift goes through 90°, i.e. the K-matrix mass. It 
does not work nearly so well using the T-matrix pole position of 1210 MeV. For Breit-Wigner resonances of constant 
width (used here), the K-matrix mass and the T-matrix pole position are the same. 

A remarkably simple pattern is apparent in Fig. 5. All states lie close to parallel straight-line trajectories. These 
trajectories extrapolate well to known states in the mass range 1200-1700 MeV. [We do not, however, attempt to 
place the rj on the 0" trajectory, since its mass is affected strongly by the instanton interaction]. The simplicity of 
Fig. 5 suggests that observed states are qq rather than hybrids. For and '5o, states expected around 1670 MeV 
are presently missing. For quantum numbers ^D^, I - \ resonances observed in the n^n"^ channel are shown; it is 
worthwhile to show these results, since the lowest two states are very well defined. 

Using PDG masses for resonances below 1.9 GeV, slopes for different quantum numbers are shown in Table 3. 
They are consistent with the same slope within statistical errors; the mean slope is 1.143 ± 0.013 GeV^ per excitation. 
The most accurate determination of the slope comes from the trajectory, which begins with the well known 
/2(1270). However, one must expect some deviations from straight lines due to local perturbations, for example (a) 
from mixing with nearby glueballs or hybrids, (b) from level repulsion between 2^ states, (c) from nearby thresholds, 
particularly in the low mass region, and (d) from variation with s of the effects of tensor and spin-orbit splitting. There 
is an 'extra' 2"+ state at 1860 MeV, which will perturb the tratectory; as discussed in Ref. [18], it is a candidate 
for a hybrid expected in that mass region. 

The 0"^ trajectory is drawn with n = 1 for /o(1370). The line does, however, go through /o(980), and it is possible 
that this state is the n = 1 ground state rather than a molecule. The error on the K-matrix mass of /o(1370) is quite 
large, so the 0"^ trajectory is defined mostly by /o(1770) and /o(2337). Because of strong mixing between scalar 
states, the interpretation of the mass range around /o(980) and /o(1370) may be complicated; it has previously been 
considered in terms of the K-matrix approach by Anisovich et al. [21]. 
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Slope 




(GeV2) 


4++ 


1.139 ±0.037 




1.107 ±0.078 




1.253 ±0.066 


'Pi 


1.131 ±0.024 


2-^ 


1.217 ±0.055 




1.120 ±0.027 




1.164 ±0.041 


0-+ 


1.183 ±0.028 


(7=1) 


1.099 ±0.029 



Table 3: Slopes of trajectories for different quantum numbers. 

In the absence of tensor and spin-orbit sphtting, ^F^, ^F^ and ^F2 states are degenerate in mass. Tensor and 
spin-orbit splitting may be assessed from the relations [22] 

AMls = [-20MC F2) - 7 mCF3) + 21 mCF4)]/ 54, (2) 
AMt = [-4M(^F2) + 7M(3F3) - 3M(^F4)]/14. (3) 

For both multiplets of Table 2, states he highest in mass, requiring significant tensor splitting. This splitting is 
in the sense predicted by one-gluon exchange. Using a linear confining potential plus one-gluon exchange, Godfrey 
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and Isgur [23] predict SMj = 9 MeV for the lower multiplet and AM^s = -20 MeV. From Table 2, one finds 
6Mt = 20 + 4.5 MeV for the lower multiplet and 7.1 + 7.3 MeV for the upper one; errors take into account observed 
correlations in fitted masses. Because we use the same centrifugal barrier for ^^"4, and states, there is almost 
no correlation of AMj or AM^ with the radius of the barrier. The observed spin-orbit splittings are SMis = 2.4 ± 4.4 
MeV for the lower multiplet and -6.3 ±5.1 MeV for the upper one; these average approximately to zero. 

In summary, the expected qq states in this mass range are observed with well determined masses and widths, 
except for the lower 0" and the upper 1"^ state, where errors are sizeable. They follow a simple pattern requiring 
approximately linear trajectories of mass squared against excitation number, with a slope of 1.143 + 0.013 GeV^ per 
excitation. In addition, there is evidence for a broad 2^ component with M - 2010 ± 25 MeV, F - 495 ± 35 MeV. 
Tensor splitting for ■^F states is, within sizeable errors, consistent with that predicted from one-gluon exchange and a 
linear confining potential. Spin-orbit sphtting of states is consistent with zero. 

1. Acknowledgement 

We thank the Crystal Barrel Collaboration for allowing use of the data. We acknowledge financial support from 
the British Particle Physics and Astronomy Research Coimcil (PPARC). We wish to thank Prof. V. V. Anisovich for 
helpful discussions. The St. Petersburg group wishes to acknowledge financial support from PPARC and INTAS grant 
RFBR 95-0267. 

References 

[1] A.V. Anisovich et al., Phys. Lett. B452 (1999) 173; Nucl. Phys. A651 (1999) 253. 

[2] A.V Anisovich et al., Phys. Lett. B468 (1999) 304 and 309; Nucl. Phys. A662 (2000) 344. 

[3] A.V Anisovich et al., Phys. Lett. B471 (1999) 271; Nucl Phys. A662 (2000) 319. 

[4] Particle Data Group (PDG), Euro. Phys. Joum. 3 (1998) 1. 

[5] D.V. Bugg, A.V. Sarantsev and B.S. Zou, Nucl. Phys. B471 (1996) 59. 

[6] A. Abele et al., Euro. Phys. J. C8 (1999) 67. 

[7] A.V. Anisovich et al., Resoiiiiiiccs in pp rjTT^Ti^n^n^ , submitted to Nucl. Phys. A. 

[8] V. V. Anisovich, D.V. Bugg, A.V. Sarantsev and B.S. Zou, Phys. Rev D50 (1994) 1972. 

[9] A.V Anisovich et al., Phys. Lett. B449 (1999) 154. 

[10] D. Barberis et al., Phys. Lett. B471 (2000) 440. 

[11] H. Baher, Analog and Digital Signal Processing, (Wiley, Chichester, 1990), p243. 

[12] D. Aide et al., Phys. Lett. B241 (1990) 600. 

[13] G.M. Beladidze et al., Zeit. Phys. C54 (1992) 367. 

[14] D.V. Bugg and B.S. Zou, Phys. Lett. B396 (1997) 295. 

[15] D.V. Bugg, Resonances around 2 GeV: qq and glueballs, Hadron99 Proceedings (to be published). 

[16] T.A. Armstrong e tal., Phys. Lett. B307 (1993) 394. 

[17] D. Barberis et al., Phys. Lett. B413 (1997) 213. 

[18] A.V Anisovich et al., Phys. Lett. B477 (2000) 19. 

[19] A.V Anisovich et al.. Data on pp ^ rfi^T^ for masses 1960 to 2410 MeV, submitted to Phys. Lett. B. 

[20] C.A. Baker et al., Phys. Lett. B467 (1999) 147. 

[21] A.V. Anisovich, V.V. Anisovich and A.V. Sarantsev, Zeit. Phys. A359 (1997) 173. 

[22] D.V. Bugg et al., J. Phys. G: Nucl. Phys 4 (1978) 1025. 

[23] S. Godfrey and N. Isgur, Phys. Rev. D32 (1985) 189. 



8 



1800 
1600 
1400 
1200 
1000 
800 
600 
400 
200 




: 900 MeV/c 












i 




I 




1 2 
M(7T7t) [ GeV ] 



1 2 
M(7T7r) [ GeV ] 



1 2 
M(7T7]) [ GeV ] 



1 2 
M(7T7;) [ GeV 1 



1 2 
M(7T7}) [ GeV ] 



1 2 
M(n7T) [ GeV ] 




Figure 1: Mass projections on to M{7rn) and M{nr]) for data at (a) and (b) 900 MeV/c, (c) and (d) 1525 MeV/c, (e) and (f) 1940 MeV/c; histograms 
show the result of the fit. Data are uncorrected for (small) variations of acceptance; these are included in the maximum likelihood fit. 
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Figure 2: Intensities fitted to individual partial waves, and corrected for all and decays. The titles in each panel indicate in descending order 
full, dashed, chain and dotted curves; o-(lOOO) denotes the broad cr (/o(400 - 1200)) in the nn amplitude. Numerical values indicate factors by 
which small intensities have been scaled so as to make them visible. 
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Figure 3: Argand diagrams for prominent partial waves; for (a) and (b), dasiied curves siiow tlie effect of removing tlie lower resonance. Crosses 
mark beam momenta of Table 1; in all cases, the amplitude moves anti-clockwise with increasing beam momentimi. 
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Figure 4: Centre of mass angular distributions for production of (a) /2(1270);;, (b) a2(1320);r, (c) ao(980);r at 1050 MeV/c; (d)-(f) show corre- 
sponding angular distributions at 1525 MeV/c. Histograms show the fit. Data and fit are uncorrected for angular acceptance of the detector; this 
acceptance is included in the maximum likelihood fit. 
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Figure 5: Plots of v. radial excitation number, n. In (b), / = 1 'D3 states are shown. Straight line trajectories are drawn in all cases with a slope 
of 1.143 GeV-^ 
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